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Abstract 

In  an  effort  to  increase  the  range  of  missiles  and  guided  bombs,  the  USAF  is 
looking  at  options  for  new  wing  configurations.  One  such  configuration  being 
considered  is  oblique  wings.  An  oblique  wing  is  a  wing  that  pivots  about  a  point  on  the 
aircraft  fuselage  thereby  having  one  side  swept  forward  and  the  other  swept  aft. 
Additional  interest  is  looking  at  a  wing  only  configuration  that  also  rotates  with  one  wing 
tip  forward  of  the  other.  Studies  have  shown  that  this  configuration  can  provide  less  drag 
for  a  given  lift  at  both  supersonic  cruise  and  subsonic  loiter  conditions.  This  experiment 
focused  on  the  low  speed  performance  of  a  missile  model  with  an  oblique  wing.  The 
wing  was  tested  at  seven  different  sweep  angles  and  at  two  different  speeds.  In  order  to 
simulate  the  missile  dropping  from  an  aircraft  the  model  was  inverted  over  a  stationary 
ground  plane  in  the  tunnel  and  tested  at  the  same  wing  obliquity  angles.  Stalling  was 
found  at  certain  conditions  including  sweep  angles  of  0,  15,  and  30  deg.  The  ground 
plane  was  shown  to  result  in  an  increase  in  lift  as  well  as  an  increase  in  drag.  The  ground 
plane  was  also  shown  to  add  more  longitudinal  stability,  thus  making  the  missile  better 
performing  when  dropped  from  an  aircraft. 
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EXPERIMENTAL  INVESTIGATION  OF  OBLIQUE  WING  AERODYNAMICS  AT 

LOW  SPEED 

I.  Introduction 

Section  1  -  Oblique  wings 

Ever  since  Chuck  Yeager  broke  the  sound  barrier  in  1947,  aircraft  designers  have 
been  looking  for  wing  designs  that  perfonn  well  at  both  subsonic  and  supersonic  speeds. 
The  general  consensus  is  that  the  only  way  to  accomplish  this  goal  was  through  wing 
morphing  of  some  kind  providing  variable  wing  geometry.  The  purpose  of  wing 
morphing  is  to  provide  a  high  aspect  ratio  for  low-speed  flight  and  for  loitering  and  also  a 
low  aspect  ratio  for  high  speed  flight. 

The  first  generation  of  variable  geometry  wings  was  the  swing  wing  designs. 
These  designs  utilized  a  pivoting  wing  that  could  swing  back  into  a  delta  configuration  at 
high  speeds.  Several  swing  wing  airplanes  were  utilized  by  the  United  States  Military  in 
the  1960s  and  1970s.  The  most  notable  and  popular  designs  were  the  F-14,  F-  111,  and 
B-l.  These  aircraft  experienced  some  success,  but  the  added  mass  and  complexity  of  the 
swing  wing  mechanism  made  their  more  widespread  use  undesirable. 

Another  variable  geometry  concept  currently  garnering  more  research  is  the 
oblique  wing  design.  The  oblique  wing  “can  vary  the  wing  sweep  with  a  single  pivot  that 
is  primarily  loaded  in  tension,  trading  aspect  ratio  for  fineness  ratio  by  sweeping  one 
wing  tip  forward  and  the  other  wing  tip  back”  [1].  Oblique  wing  aircraft  have  fewer 
moving  parts  than  swing  wing  aircraft,  thus  making  them  lighter  and  less  complex. 

The  first  known  oblique  wing  design  was  the  Blohm  and  Voss  P-202,  designed  in 
Gennany  by  Richard  Vogt  in  1942.  The  P-202  can  be  seen  in  Fig.  1.  The  aircraft  was  a 
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single  engine  fighter  that  featured  a  top  mounted  wing  designed  to  rotate  in  order  to  slow 
the  onset  of  Mach  divergence  [1].  In  1952,  R.T.  Jones  proved  that  for  any  flight  Mach 
number  the  minimum  drag  for  a  given  lift  could  be  achieved  by  an  oblique  swept  wing 
with  an  elliptic  planform  [2],  Other  wind  tunnel  tests  have  shown  that  oblique  wings  are 
very  effective  at  reducing  wave  drag  at  supersonic  speeds  [3],  More  recently  NASA  test 
flew  a  full  scale  demonstrator,  the  AD-1.  The  AD-1  flew  for  79  hours  and  performed 
well,  although  there  were  serious  handling  qualities  issues  at  wing  sweeps  of  greater  than 
45  deg  [4].  Fig.  2  shows  the  AD-1. 


’  Hd-V  l\ 202 


Figure  1.  Blohm  and  Voss  P-202  [1], 


Figure  2.  NASA  AD-1  [4]. 
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Section  2  -  Ground  effects 

Ever  since  the  advent  of  powered  flight,  pilots  have  experienced  a  “pillow  effect” 
when  flying  close  to  the  ground.  This  phenomenon  is  most  felt  during  the  most  critical 
and  dangerous  phases  of  flight,  takeoff  and  landing.  Most  aircraft  experience  a 
significant  increase  in  lift  and  a  decrease  in  induced  drag.  Pilots  must  manually  account 
for  this  phenomenon  because  aircraft  are  not  designed  to  deal  with  this  condition  [5]. 

An  aircraft  is  normally  considered  to  be  in-ground  effect  (IGE)  when  it  is  flying 
within  one  wingspan,  b,  of  the  ground.  When  the  aircraft  is  IGE  the  induced  drag  is 
usually  reduced  because  of  the  interaction  between  the  wingtip  vortices  and  the  ground 
[6].  While  the  aircraft  is  flying  out-of-ground  effect  (OGE)  the  vortices  are  shed  off  the 
wing  in  a  cylindrical  pattern.  During  IGE  flight  the  vortices  begin  to  flatten  out  because 
the  presence  of  the  ground  prevents  them  from  forming  cylindrically.  This  causes  an 
increase  in  effective  aspect  ratio  and  wingspan.  Since  aspect  ratio  and  induced  drag  are 
inversely  proportional,  an  aircraft  flying  IGE  will  experience  a  lower  induced  drag  and 
therefore  a  decreased  total  drag,  assuming  profile  drag  remains  constant  [6], 

Along  with  the  reduction  in  induced  drag,  aircraft  flying  IGE  can  be  expected  to 
have  a  higher  lift  and  pitching  moment.  An  important  measure  of  aircraft  aerodynamic 
efficiency  is  the  lift  to  drag  ratio,  L/D.  Since  an  aircraft  IGE  has  a  higher  lift  and  lower 
drag  it  performs  much  more  efficiently  than  an  identical  one  flying  OGE  at  the  same 
flight  conditions.  This  fact  led  Russian  designer  Rostislav  Alexiev  to  develop  what  he 
called  an  ekranoplan.  He  designed  the  KM  Caspian  Sea  Monster,  a  550  ton  vehicle  for 
fast  transportation  and  heavy  cargo  carrying  over  water.  The  Soviets  designed  several 
more  ekranoplans  before  the  program  was  cancelled  due  to  lack  of  funding  in  the  late 
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1980s  [7].  Boeing  is  also  currently  examining  the  viability  of  a  similar  vehicle  to  the 
Russian’s  designs  capable  of  carrying  an  entire  Anny  Division  or  up  to  17  M-l  tanks  [8]. 
Section  3  -  Research  objectives  and  approach 

The  goals  of  this  research  were  to  determine  the  aerodynamic  properties  of  a 
missile  with  an  oblique  wing  installed  and  to  see  how  the  model  performed  in  ground 
effect.  The  AFIT  low  speed  wind  tunnel  was  used,  both  without  a  ground  plane  for  OGE 
tests  and  with  the  ground  plane  installed  to  determine  IGE  behavior.  The  IGE  tests  were 
run  with  the  missile  inverted  above  the  ground  plane,  with  the  ground  plane  simulating  an 
aircraft.  The  IGE  tests  were  accomplished  to  detennine  any  undesirable  effects  of  an 
oblique  winged  weapon  shortly  after  separation  from  an  aircraft.  The  tests  were 
performed  at  low  speed  for  applicability  to  low  speed  guided  bomb  drops 
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II.  Literature  Review 


Section  1:  Oblique  Wings 

Oblique  wings  have  been  studied  since  World  War  II  due  to  their  promise  of 
aerodynamic  efficiency  for  aircraft  flying  at  both  transonic  and  supersonic  speeds  [1].  In 
1958,  R.  T.  Jones  noted  that  wave  drag  and  induced  drag  could  be  minimized  by  a 
variable  sweep  oblique  wing  with  an  elliptical  lift  distribution  [2].  At  low  speed  the  most 
efficient  wing  has  a  high  aspect  ratio  and  is  elliptically  loaded  due  to  the  fact  that  induced 
drag  is  inversely  proportional  to  aspect  ratio.  At  supersonic  speeds  drag  on  an  aircraft  is 
dominated  by  wave  drag.  One  major  advantage  of  oblique  wings  for  supersonic  flight  is 
that  for  equivalent  span,  sweep,  and  volume  they  distribute  the  lift  over  twice  the  length 
of  a  more  conventional  swept  wing  planform.  Fig.  3  illustrates  this  point.  This  reduces 
lift  dependent  wave  drag  by  a  factor  of  four  and  volume  dependent  wave  drag  by  a  factor 
of  16  [1], 


Figure  3.  Illustration  of  2X  length  for  oblique  wings  [1], 
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Jones  also  showed  that  the  induced  drag  of  an  oblique  wing  at  optimal  sweep  is 
half  of  that  for  a  delta  wing  of  the  same  span  [2],  Fig.  4  shows  the  induced  drag  of  an 
oblique  wing  and  delta  wing. 


Figure  4.  Comparison  of  induced  drag  for  oblique  and  delta  wing  [9] 

One  of  the  main  problems  associated  with  the  use  of  oblique  wings  is 
controllability.  At  low  sweep  angles,  an  oblique  winged  aircraft  can  be  controlled  much 
like  a  conventionally  winged  vehicle.  In  fact,  using  the  NASA  Langley  tunnel,  Campbell 
noted  that  at  sweep  angles  of  40  deg  or  less  there  were  no  serious  controllability  issues 
encountered  [10].  However,  he  also  noted  that  at  sweep  angles  above  60  deg  ailerons 
become  unsatisfactorily  weak.  He  theorized  that,  “The  aileron  rolling  effectiveness  was 
not  reduced  by  skewing  the  wing  from  0  deg  to  40  deg  because  the  damping  in  roll 
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decreased  approximately  the  same  amount  as  the  aileron  rolling  moments  [10].”  He  also 
noted  that  at  a  sweep  angle  of  40  deg  a  large  pitching  moment  was  present,  although  this 
was  not  experienced  during  flight  testing.  This  is  because  “the  lift  forces  on  the  wing 
produced  by  rolling  introduced  pitch  moments  were  equal  to  and  opposite  of  the  aileron 
pitching  moments  [10].” 

Kempel  et  al.  [11]  performed  flying  qualities  testing  on  the  NASA  F-8  oblique 
wing  research  aircraft.  They  noted  that  in  general,  as  the  sweep  angle  and  dynamic 
pressure  increased,  pilot  ratings  decreased.  They  also  noticed  that  in  open  loop 
configurations  there  was  significant  pitch-to-roll  coupling  and  an  unacceptable  amount  of 
pitch-to-sideforce  coupling.  At  high  dynamic  pressures  one  pilot  even  described  the 
lateral  response  to  a  pitch  up  command  as  “scary”  [11].  This  clearly  means  that 
sophisticated  multi-input-multi-output  control  systems  will  need  to  be  developed  for 
oblique  wing  aircraft  before  they  begin  flying  on  a  larger  scale. 

Section  2:  Ground  Effects 

Ground  effect  is  a  phenomenon  caused  by  the  presence  of  a  boundary  near  a  wing 
altering  the  airflow  around  a  wing.  This  altered  airflow  causes  an  increase  in  the  lift  and 
a  decrease  in  the  induced  drag  of  a  wing.  The  closer  the  wing  is  to  the  boundary,  the 
more  exaggerated  this  effect  becomes.  This  is  caused  by  the  air  flowing  underneath  the 
wing  not  being  able  to  expand  as  it  would  in  free  air  [12].  The  difference  in  the  airflow 
around  a  wing  OGE  and  one  flying  IGE  can  be  seen  in  Fig.  5.  Halloran  and  O’Meara 
[12]  describe  the  effect  in  terms  of  pressure: 

“In  terms  of  the  total  pressure  of  the  flow,  the  additional  lift  is  due  to  a  rise  in 
static  pressure  under  the  wing.  The  total  pressure  of  the  flow  field  can  be  divided  between 
the  static  pressure  (surface  pressure)  and  dynamic  pressure  (the  pressure  associated  with 
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velocity).  As  the  total  pressure  remains  constant  throughout  the  flow  field,  the  sum  of  the 
static  and  dynamic  pressure  must  also  remain  constant.  As  the  flow  is  forced  into  the 
region  between  the  wing  and  the  boundary,  the  decrease  in  dynamic  pressure  is 
transformed  into  a  rise  in  the  static  pressure.  This  rise  in  the  static  pressure  is  often 
referred  to  as  ‘ram  pressure’.  The  resulting  altered  pressure  distribution  causes  a  net 
increase  in  the  lift  and  a  change  to  many  of  the  other  aerodynamic  characteristics  of  the 
wing  [12].” 


Figure  5.  Airflow  around  wing  OGE  and  wing  IGE  [12], 


Also,  there  is  not  enough  space  for  the  vortices  to  fully  develop  when  a  wing  is 
approaching  the  ground.  Therefore  the  amount  of  leakage  of  pressure  from  the  lower 
side  to  the  upper  side  is  less  and  the  vortices  become  weaker.  The  vortices  are  also 
pushed  outward  by  the  ground,  causing  the  effective  aspect  ratio  to  become  higher  than 
the  geometric  aspect  ratio  [12], 

The  “pillow  effect”  caused  by  the  increase  in  lift  an  aircraft  experiences  when 
close  to  the  ground  has  been  experienced  by  pilots  since  the  beginning  of  manned  flight. 
In  the  early  1920s,  a  German  engineer  named  Carl  Wieselsberger  used  Prandl’s  three- 
dimensional  lifting  line  theory  to  develop  an  equation  for  the  reduction  in  induced  drag 
for  an  aircraft  flying  near  the  ground  [13].  This  equation  has  been  used  as  the  standard 
and  has  been  verified  by  flight  tests  throughout  the  years  [5].  The  Wieselsberger 
equation  relates  the  induced  drag,  C’d,,  to  the  h/b  of  the  wing.  Fig.  6  shows  a  plot  of 
Wieselsberger’ s  equation,  with  Coi  normalized  to  the  value  OGE. 
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Mieselsberger 


Figure  6.  Graph  of  Weiselsberger’s  induced  drag  formula  [7], 

McCormick  took  a  different  approach  when  developing  equations  for  lift  and  drag 
to  determine  ground  roll  of  an  aircraft.  He  used  the  lifting  line  theory  and  the  Biot-Savart 
law  to  come  up  with  an  equation  for  the  ratio  of  C’d,  IGE  to  Cdi  OGE.  He  modeled  the 
wing  as  a  horseshoe  vortex  with  the  ground  simulated  by  having  the  vertical  velocities 
canceling  each  other  out.  He  used  the  height  as  the  difference  between  the  reflection 
plane  and  the  vortices.  Next,  the  Biot-Savart  law  estimated  the  velocities  induced  at  a 
certain  point  from  every  horseshoe  vortex.  This  allowed  him  to  come  up  with  the 
following  relationship  [14]: 


CnXiGE) 

CDi{OGE) 


K  2 


16(-) 

b 


K  2 


1+(h} 
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Fig.  7  shows  a  plot  of  McCormick’s  induced  drag  factor.  It  can  be  seen  that 
although  McCormick  and  Wieselsberger  took  different  approaches  they  both  came  up 
with  very  similar  results. 
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Figure  7.  McCormick’s  induced  drag  factor  [14] 

Some  of  the  first  ground  effects  tests  were  accomplished  by  Raymond  at  the 
Massachusetts  Institute  of  Technology  in  1921  [15].  He  installed  a  flat  plate  in  the  tunnel 
and  tested  three  different  airfoils  at  varying  heights  above  the  plate.  He  also  attempted  to 
use  reflection  to  create  an  imaginary  ground  plane.  Both  sets  of  tests  showed  that  near 
the  ground  a  wing  will  experience  an  increase  in  lift  and  a  decrease  in  induced  drag  [15]. 

Raymond’s  technique  of  using  a  fixed  ground  plane  height  and  model  has 
become  known  as  static  wind  tunnel  testing.  Different  ground  heights  were  tested  by 
moving  the  ground  plane  either  closer  or  farther  from  the  model.  Pilots  flew  ground 
effect  tests  known  as  ‘fly-by’  patterns  in  an  attempt  to  validate  the  tunnel  results.  Angle 
of  attack  and  altitude  were  kept  constant  to  find  the  location  of  the  ground  effect  region 
[5].  This  technique  fails  to  capture  the  true  dynamics  of  an  aircraft  as  it  approaches  a 
runway  because  the  aircraft  has  no  sink  rate.  Schweikhard  developed  a  method  in  1967 
to  measure  the  ground  effect  of  an  aircraft  in  flight  [16].  Test  pilots  would  fly  the  aircraft 
at  a  set  angle  of  attack  and  throttle  setting,  but  would  allow  the  sink  rate  to  vary.  This 
technique  ensured  that  lift,  drag,  and  pitch  moment  were  constant  as  the  ground 
approached.  Engineers  then  measured  any  change  in  flight  path  angle,  velocity,  or 
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control  surface  deflection  as  the  aircraft  entered  the  ground  effect  region.  The  engineers 
discovered  this  technique  saved  time  and  required  less  data  analysis  than  previous 
techniques  [16]. 

Since  wind  tunnel  testing  saves  a  lot  of  time  and  money  over  flight  testing, 
engineers  developed  new  ways  to  dynamically  test  ground  effects  in  the  wind  tunnel.  A 
dynamic  wind  tunnel  test  moves  the  model  mechanically  or  manually  towards  the  ground 
plane  to  better  simulate  take  off  or  landing  approach  [5],  Chang  et  al.  noted  the  need  for 
dynamic  testing  because  of  the  disparity  between  static  tests  and  flight  test  data  [17].  He 
tested  several  wings  and  aircraft  models  both  statically  and  dynamically  in  the  wind 
tunnel.  He  and  Baker  et  al.  noted  that  at  h/b  <  0.4,  the  static  wind  tunnel  results  over 
estimated  the  change  in  lift  due  to  ground  effect  [18], 

Most  of  the  time  an  aircraft  flying  in  the  ground  effect  region  will  experience  an 
increase  in  lift  and  a  decrease  in  drag,  but  this  is  not  always  the  case.  Lee  et  al.  [19] 
experienced  an  increase  in  lift  and  drag  as  height  above  the  ground  plane  decreased 
during  their  tests  of  a  60  deg  delta  wing,  FI 06,  and  XB-70-1  models.  They  tested  h/b 
from  0.2  to  1.6  and  Reynolds  numbers  from  3x  1  O'  to  7.5xl05  for  all  three  models.  Jones 
[5]  extrapolated  the  data  from  the  Cd  vs  h/b  and  Cl  vs  h/b  to  analyze  trends  in  L/D.  He 
noted  the  60  deg  delta  wing  experienced  a  small  decrease  in  L/D  and  that  the  F-106  and 
XB-70-1  both  had  a  decrease  and  a  small  increase  in  L/D  depending  on  h/b.  The  XB-70- 
1  experienced  a  downward  trend  in  CD  between  h/b  and  0.3  and  0.2,  which  explains  the 
increase  in  L/D  [5]. 

Curry  [20]  also  tested  the  F-16XL  for  ground  effects  using  both  a  wind  tunnel  and 
flight  testing.  He  noted  that  CD  increased  as  h/b  decreased.  He  stated  that  for  an  aircraft 
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flying  close  to  the  ground,  an  increase  or  decrease  in  drag  is  possible.  An  increase  in 
drag  was  also  noted  on  the  Tu-144  supersonic  transporter  when  flying  close  to  the  ground 
by  Curry  and  Owens  [21]. 

One  main  concern  when  using  a  ground  plane  in  a  wind  tunnel  is  the  boundary 
layer  build-up  on  the  top  surface  of  the  ground  plane.  Any  surface  that  has  direct  contact 
with  a  moving  fluid  will  experience  a  boundary  layer  forming.  This  boundary  layer  can 
cause  inaccurate  data  when  a  ground  plane  is  installed.  The  solution  to  this  problem  is  a 
boundary  layer  removal  system  [5]. 

The  most  common  way  to  remove  a  boundary  layer  is  through  the  use  of  a 
moving  belt  system.  The  belt  moves  in  the  tunnel  at  the  same  speed  as  the  airflow  thus 
better  simulating  an  aircraft  flying  over  a  stationary  ground.  However,  several  studies 
have  shown  that  a  moving  belt  is  not  required  to  get  accurate  results  at  all  flight 
conditions.  It  has  been  shown  that  the  need  for  a  moving  belt  is  dependent  on  the 
maximum  CL  of  the  aircraft.  Turner  [22]  determined  the  conditions  under  which  an 
endless  belt  system  would  be  needed  for  accurate  results.  He  showed  that  the  necessity 
of  a  moving  belt  depended  on  the  spanwise  lift  coefficient  and  height  above  the  ground 
[22],  Fig.  8  shows  the  criteria  Turner  established  for  when  an  endless  belt  ground  plane 
was  needed. 
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Figure  8.  Conditions  needing  an  endless  belt  ground  plane  [22] 

The  shaded  region  represents  the  tunnel  conditions  tested  in  the  current 
study.  Therefore,  based  on  Turner’s  study  [22],  a  moving  belt  was  not  required  for  this 
study,  and  one  was  not  used. 

Kemmerly  and  Paulson  [23]  did  another  study  on  the  necessity  of  the  use  of  a 
moving  plane.  They  found  that  if  equation  2  was  met,  then  the  use  of  a  moving  ground 
plane  is  needed. 


h/b 

"c7 


<0.05 


(2) 


Table  1  shows  the  values  used  in  this  study.  It  can  be  seen  that  none  of  the 
configurations  gives  a  value  of  less  than  0.05,  so  once  again  it  is  confirmed  that  a  moving 
ground  plane  is  not  needed  for  boundary  layer  removal.  You  can  also  see  that  if  equation 
(2)  was  plotted  similarly  to  Fig.  8,  it  would  give  almost  equal  results. 
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Table  1.  Justification  for  use  of  stationary  ground  plane 


Sweep  angle,  X 
(deg) 

Cl  max 

h/b  min 

(h/b)/CLmax 

<0.05 

0 

0.78 

0.071 

0.09 

No 

15 

0.3 

0.074 

0.25 

No 

30 

0.45 

0.082 

0.18 

No 

45 

0.7 

0.101 

0.14 

No 

60 

0.55 

0.143 

0.26 

No 

75 

0.3 

0.276 

0.92 

No 

90 

0.22 

0.625 

2.84 

No 
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III.  EXPERIMENTAL  SETUP 


Section  1  -  Wind  Tunnel 

The  AFIT  3  ’  X  3  ’  wind  tunnel  was  used  for  this  experiment.  The  tunnel  was 
manufactured  by  New  York  Blower  Company.  The  wind  tunnel  is  powered  by  an 
ACF/PLR  class  IV  fan  with  a  Toshiba  Premium  Efficiency  (EQP  III)  fan  motor.  Both  the 
fan  and  the  motor  are  controlled  by  a  Siemens  (13710)  Adjustable  Frequency  Tunnel 
Controller.  The  tunnel  has  an  operating  speed  of  1785  rpm  and  has  a  maximum  speed  of 
150  mph.  Table  2  shows  the  fan  and  motor  specifications: 


Table  2.  Fan  and  Motor  Specifications  [24]. 


Motor 

Controller 

3  phase  induction 

1785  RPM  operating  speed 

Maximum  theoretical  speed  -  150  mph 

Maximum  tested  speed  -  148  mph 

200  brake  horsepower 

250  maximum  horsepower 

230/460  volts 

460  volts 

444/222  amps 

315  amps 

60  Hz 

4  poles 

The  tunnel  is  an  Eiffel-type  with  an  open  circuit  configuration  and  a  closed  test 
section  44”  X  3 1”  in  size.  The  fan  is  at  the  end  of  the  tunnel  and  pulls  in  ambient  air 
through  a  122”  W  X  1 1 1”  H  X  70”  L  intake  plenum.  The  plenum  consists  of  four  steel 
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mesh  anti-turbulence  screens  with  %”  aluminum  honeycomb  flow  straighteners  with  a 
minimum  aspect  ration  of  15  to  give  the  tunnel  good,  laminar  streamlines  [25].  The  flow 
travels  to  the  test  section  through  a  95.5”  long  converging  section  with  a  contraction  ratio 
of  9.5:1.  The  height  of  the  tunnel  is  111”  through  the  anti-turbulence  screens  and 
contracts  to  the  test  section  height  of  31”.  Fig.  9  shows  the  dimensions  of  the  converging 
section  of  the  tunnel  and  Fig.  10  displays  a  schematic  of  the  tunnel. 


Figure  9.  Wind  tunnel  converging  section  dimensions. 


Eifcaurt 


Figure  10.  AFIT  wind  tunnel  schematic. 
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The  test  section  has  an  octagonal  shape  to  reduce  corner  interference  [24].  It  has 
dimensions  of  3 1”H  X  44”W  X  72”L.  The  test  section  is  shown  in  Fig.  1 1 .  The  test 
section  sting  mechanism  is  remotely  controlled  and  can  vary  angles  of  attack  from  -25  to 
25  degrees. 


Wind 

Input 


Wind 

Exhaust 


Figure  11.  Wind  tunnel  test  section. 


The  tunnel  diverges  after  the  test  section  and  the  air  exits  through  a  vertical 
exhaust  plume.  There  is  a  fence  at  the  end  of  the  tunnel  to  prevent  any  debris  from 
impacting  the  fan  or  motor. 

Section  2  -  Balance 

The  balance  used  for  this  experiment  was  AFIT’s  eight  component  25-lb  balance 
manufactured  by  the  Able  Corporation  (0.5  MK  VI).  It  is  an  internal  balance  and  uses 
strain  gauge  rosettes  to  measure  loads  in  the  various  directions.  The  maximum  loads  of 
this  balance  in  the  each  specific  direction  are  shown  in  Table  3. 
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Table  3.  Maximum  loads  for  251b  balance. 


Component 

Maximum 

Normal  Force  (Nl) 

25  lbs 

Pitch  Moment  (N2) 

25  in-lbs 

Side  Force  (SI) 

15  lbs 

Yaw  Moment  (S2) 

15  in- lbs 

Axial  Force  (Al) 

30  lbs 

Roll  Moment  (LI) 

15  in-lbs 

The  balance  measures  forces  in  terms  of  changes  in  voltage.  LabView,  which  is 
used  to  write  data  acquisition  and  calibration  programs,  collects  the  data  and  converts  the 
voltages  to  pounds  force  for  use  in  calculating,  lift,  drag,  moments,  and  other  parameters. 
Section  3  -  Missile  Model  and  Wing 

An  aluminum  missile  model  28.44  in  long  with  a  projected  diameter  of  2  in  was 
used  in  this  experiment.  The  missile  has  four  identical  tail  (ins,  two  horizontal  and  two 
vertical,  made  from  a  symmetric  airfoil.  The  bare  missile  weighs  7.33  lbf.  The  missile  is 
shown  in  Fig.  12. 


< - 

Airflow  direction 

◄ - 

◄ - 


Figure  12.  Missile  model  mounted  in  tunnel. 
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The  wing  section  used  for  the  wing  is  shown  in  Fig.  13.  It  has  a  chord  of  2  in,  a 
span  at  A=0  of  17.5  in,  has  no  twist  or  dihederal,  but  does  have  a  positive  camber.  This 
airfoil  was  chosen  in  order  to  be  able  to  compare  the  results  of  this  experiment  to 
previous  data  using  this  missile  model.  Corneille  [26]  tested  a  joined  wing  configuration 
to  determine  if  it  was  beneficial  over  a  missile  with  a  single  wing.  Dike  [25]  then  tested 
the  joined  wing  as  it  morphed  out  to  full  extension.  In  both  cases  the  missile  model  was 
the  same  one  used  in  this  study  and  the  joined  wings  tested  had  the  same  basic  airfoil 
shape.  For  this  experiment  the  chord  was  doubled  in  order  to  add  more  stiffness  and 
strength  to  support  the  loading  of  the  entire  span. 


H -  2  in  - ►! 


Figure  13.  Airfoil  profile. 

In  this  study,  the  sweep  angle  X  is  measured  with  the  wing  perpendicular  to  the 
missile  being  0  deg  and  the  wing  parallel  to  the  missile  being  90  deg.  The  convention 
used  to  define  the  sweep  angle  is  shown  in  Fig.  14. 
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Figure  14.  Sweep  angle  convention  used 

The  wing  was  first  drawn  in  SolidWorks  and  then  sent  to  AFIT’s  3-D  rapid 
prototyping  machine.  The  Stratasys  EDEN  333  rapid  prototyping  machine  lays  down 
0.0006  in  layers  of  LTV  plastic  (also  called  photopolymer  plastic)  and  a  gel-like  UV 
plastic  for  support  materials  using  eight  small  jets.  The  eight  jets  traverse  the  printed 
region  in  2  in  strips  and  then  a  UV  light  cures  the  plastics  simultaneously  [27].  The  Full 
Cure  700  series  photopolymer  plastic  model  material  can  be  machined,  drilled,  and 
painted  [28].  The  SolidWorks  drawing  of  the  wing  can  be  seen  in  Fig.  15. 
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Figure  15.  SolidWorks  drawing  of  wing. 

The  wing  was  installed  on  the  missile  at  two  different  locations.  The  tests  were 
run  with  the  leading  edge  of  the  wing  installed  4  in  from  the  nose  of  the  missile,  referred 
to  as  front  wing  in  this  study,  and  were  also  run  with  the  wing  at  8  14  in  from  the  nose, 
referred  to  as  middle  wing,  to  try  and  determine  if  the  wing  was  getting  airflow  disturbed 
by  the  nose  of  the  missile. 

Section  4  -  Ground  Plane 

A  ground  plane  was  installed  in  the  tunnel  to  simulate  the  missile  model 
separating  from  the  aircraft.  The  ground  plane  was  made  of  two  steel  plates  and  eight 
cylindrical  steel  legs.  The  dimensions  are  shown  in  Table  4,  and  a  picture  of  the  ground 
plane  is  shown  in  Fig.  16. 
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Figure  16.  Ground  plane  [5]. 


Table  4.  Ground  plane  dimensions  [5]. 


Plate 

Thickness,  in 

0.25 

Diameter/width,  in 

35.31 

Max  length,  in 

44.31 

Legs 

Diameter 

1.5 

Length,  in 

Height  1 

9.77 

Height  2 

12.17 

Height  3 

12.97 

Height  4 

13.77 
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The  circular  plate  is  the  same  size  as  the  bottom  plate  of  the  test  section,  which  is 
rotated  to  simulate  yaw  angles.  Mounting  the  circular  plate  in  this  manner  allows  models 
being  tested  with  the  ground  plane  to  be  run  at  the  same  yaw  angles  as  those  tested 
without  the  ground  plane.  There  is  also  a  slit  in  the  circular  plate  to  reduce  interference 
caused  by  the  sting  mechanism  rotating  to  simulate  angle  of  attack.  The  front  piece  is 
rounded  and  beveled  and  provides  a  straight  leading  [5],  Fig.  17  shows  the  two  ground 
plane  plates  separated  and  Fig.  18  shows  the  ground  plane  leading  edge. 


Figure  17.  Top  view  of  ground  plane  with  pieces  separated  [5]. 
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Figure  18.  Ground  plane  leading  edge  [5]. 


Section  5  -  Experimental  Procedure 

The  first  step  in  taking  data  was  to  perform  a  balance  calibration  using  static 
weights.  Known  loads  were  applied  in  the  different  directions  to  the  balance  and  the 
calibration  constants  were  calculated  by  the  software  by  matching  the  loads  on  the 
balance  to  the  loads  recorded  in  the  software.  It  was  ensured  that  the  registered  voltages 
corresponded  linearly  to  the  increases  of  the  applied  weights  so  that  linearity  could  be 
verified.  Lab  View  Virtual  Instrument  was  used  to  control  all  tunnel  parameters.  These 
parameters  include  angle  of  attack,  yaw  angle,  and  tunnel  velocity.  In  order  to  provide 
analog  backups  for  the  angle  of  attack  and  sideslip  angle,  there  are  also  optical  encoders 
mounted  on  the  sting.  For  the  velocity,  a  pressure  transducer  and  pitot-static  tube  were 
used  as  the  analog  measurement.  These  measurements  were  used  as  the  main  guide  for 
tunnel  velocities  throughout  all  of  the  testing  [5]. 
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The  balance  data  were  stored  as  two  normal  force  components  (Ni  and  N?),  two 
side  force  components  (Si  and  S2),  an  axial  force  component  (Ai),  and  roll  moment  (£1). 
The  strain  gage  rosette  experienced  a  continuous  flow  of  voltage,  with  the  resistance 
being  constantly  measured.  Any  increase  in  the  length  of  the  wire  was  caused  by  an 
increase  in  the  resistance.  The  output  voltages  caused  by  the  increased  resistance  were 
equated  to  strain  and  ultimately  force  through  a  series  of  calibration  equations.  A 
conventional  coordinate  system  was  used  with  the  positive  x-direction  pointing  towards 
the  intake,  the  positive  y-direction  pointing  right  when  looking  in  the  positive  x-direction, 
and  the  positive  z-direction  pointing  down  [5].  Fig.  19  is  an  illustration  of  the  coordinate 
system  used. 
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Figure  19.  Tunnel  coordinate  system  [29]. 

The  missile  and  wing  were  tested  in  two  different  flight  conditions.  The  first  was 
with  the  wing  mounted  on  top  of  the  missile  to  measure  forces  out-of-ground-effect 
(OGE).  The  proposed  test  conditions  include  runs  at  speeds  of  100,  130,  and  145  mph. 
At  each  speed  the  sweep  angle,  X,  was  varied  from  0  to  90  deg,  in  15  deg  increments 
while  angle  of  attack,  a,  was  varied  from  -4  to  15  degrees  in  2  deg  increments  until  the 
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final  1  deg  increment  from  14  deg  to  15  deg.  On  the  first  run,  with  k  =  0  deg  and  speed 
of  100  mph,  the  model  began  to  oscillate  violently  up  and  down  at  a  =  4  deg  so  the  run 
was  aborted.  The  decision  was  then  made  to  run  the  tests  at  50  and  60  mph,  with  all 
other  parameters  as  proposed.  A  tare  run  was  also  completed  to  determine  the  effects  of 
the  missile’s  and  wing’s  static  weight  on  the  balance. 

Once  the  OGE  tests  were  complete,  the  missile  was  rotated  180  degrees  and  the 
wing  mounted  inverted  underneath  the  missile  and  above  the  ground  plane.  This  was 
done  to  simulate  launch  from  an  aircraft,  and  to  detennine  the  forces  caused  by  in- 
ground-effects  (IGE).  The  original  plan  was  to  run  the  tests  with  all  four  ground  plane 
heights.  After  running  the  tests  at  the  first  two  heights  it  was  determined  that  testing  at 
the  other  two  ground  plane  heights  would  not  be  possible  due  to  inadequate  clearance. 
With  the  first  ground  plane,  a  was  varied  from  -4  to  15  deg  to  match  the  testing  without 
the  ground  plane.  With  the  second  plane  a  could  only  be  varied  from  -2  to  5  deg  due  to 
limited  clearance  and  the  possibility  of  interference  between  the  model  and  the  ground 
plane.  Table  5  shows  the  experimental  test  matrix. 


Table  5.  Experimental  test  matrix. 


Tunnel  Speed 

OGE 

IGE  Plane  1 

IGE  Plane  2 

(mph) 

(deg) 

(deg) 

(deg) 

50 

-4<  a<15 

-4<  a<15 

-2<  a<5 

60 

-4<  a<15 

-4<  a<15 

-2<  a<5 

A  very  important  parameter  when  calculating  IGE  forces  and  moments  is  h/b, 
with  h  being  the  height  above  the  ground  plane  measured  from  the  plane  to  the  nearest 
point  on  the  missile  at  0  deg  angle  of  attack,  and  b  being  the  wing  span.  When  dealing 
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with  oblique  wings  h/b  is  not  constant,  because  as  A  varies  so  does  b.  As  the  wing  is 
swept  to  higher  sweep  angles,  the  wing  span  decreases  to  a  minimum  of  2-in,  the  chord, 
when  the  wing  is  at  a  90  deg  sweep  angle.  Table  6  shows  the  relationships  between  A  and 
h/b. 


Table  6.  h/b  for  varying  A. 


h/b 

A (deg) 

b  (in) 

Plane  1 

Plane  2 

0 

17.5 

0.221 

0.071 

15 

16.9 

0.229 

0.074 

30 

15.16 

0.248 

0.082 

45 

12.37 

0.313 

0.101 

60 

8.75 

0.443 

0.143 

75 

4.53 

0.855 

0.276 

90 

2 

1.94 

0.625 

During  the  OGE  testing  the  normal  convention  for  coordinate  systems  was  used. 
Lift  was  in  the  negative  Z-direction,  drag  in  the  negative  X-direction,  and  a  positive 
pitching  moment,  measured  about  the  center  of  gravity,  causing  the  nose  of  the  missile  to 
move  in  the  negative  Z-direction.  A  positive  angle  of  attack  also  causes  the  nose  to  move 
in  the  negative  Z-direction.  Fig.  20  shows  this  configuration. 
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Figure  20.  Force  directions  for  OGE  testing  at  k  =  0. 

When  the  model  was  tested  with  the  ground  plane,  lift  and  pitching  moment  now 
acted  in  the  opposite  direction.  Now  the  lift  the  wing  created  acted  in  the  positive  Z- 
direction.  Also,  since  lift  and  positive  angle  of  attack  should  act  in  the  same  direction, 
angle  of  attack  now  caused  the  nose  to  move  in  the  positive  Z-direction.  The  drag  acts  in 
the  same  direction  for  both  configurations.  Fig.  21  shows  the  direction  in  which  the  IGE 
forces  acted. 
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Figure  21.  Force  directions  for  IGE  testing  at  k  =  0. 

To  account  for  the  directions  the  forces  acted  in,  some  of  the  data  had  to  be 
altered  before  they  could  be  analyzed.  For  the  analysis,  the  lift,  angle  of  attack,  and  pitch 
moment  that  was  registered  by  the  balance  was  multiplied  by  negative  one.  This  ensured 
that  all  of  the  forces  and  moments  measured  in  the  OGE  and  IGE  tests  were  acting  in  the 
same  direction  so  they  could  be  compared. 
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Section  6  -  Data  Analysis 

A  data  acquisition  program  was  set  up  within  the  control  computer  to  store  the 
data  in  a  tab  delimited  text  file  at  a  sampling  rate  of  2  Hz.  First,  the  flow  velocity  was 
increased  from  zero  up  to  the  desired  velocity.  After  ensuring  the  balance  was  accurately 
taking  data,  the  model  was  dropped  to  its  lowest  a  value  and  data  was  acquired  for  ten 
seconds.  The  angle  of  attack  was  then  increased  by  2  deg  and  held  for  an  additional  ten 
seconds.  This  was  repeated  until  an  angle  of  attack  of  14  deg  was  reached,  at  which  time 
it  was  increased  by  one  degree  to  15  degrees.  For  the  tests  with  the  second  ground  plane 
a  was  increased  in  one  degree  increments  until  live  degrees  was  met. 

A  MATLAB'  code  written  by  Deluca  [24]  and  Gebbie  [30]  and  altered  by  Jones 
[5]  and  In  [29]  was  used  to  reduce  the  acquired  force  and  moment  data.  The  reduction 
program  first  loads  the  appropriate  tare  and  experiment  text  files.  It  then  combines  the 
readings  for  each  a  and  averages  them  into  a  single  test  point  for  each  a.  Before  the  data 
were  exported  as  aerodynamic  coefficients,  the  testing  conditions,  balance  interactions, 
and  blockage  corrections  were  calculated  [5]. 

After  the  MATLAB  code  reduced  the  data,  an  output  file  was  created  in 
EXCEL  format  that  consisted  of  Mach  number,  Reynolds  number,  dynamic  pressure, 
velocity,  angle  of  attack,  lift,  drag,  roll  moment,  pitch  moment,  yaw  moment,  and  side 
force  for  every  angle  of  attack  tested  [5],  See  Appendix  A  for  a  sample  calculation  of  the 
data  reduction. 
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IV.  RESULTS  &  ANALYSIS 


Verification  of  balance  and  data  reduction  program 

Several  tests  were  run  to  verify  that  the  balance  was  properly  calibrated  and 
recording  accurate  raw  data  and  to  ensure  the  MATLAB  program  was  correctly 
reducing  the  raw  data.  The  missile  model  was  run  without  the  wing  installed  at  identical 
conditions  to  the  ones  run  by  Corneille  [26]  and  Dike  [25].  The  lambda  UCAV  model 
was  also  tested  at  the  same  conditions  run  by  In  [29].  Tables  7  and  8  show  a  summary  of 
flight  conditions  tested  using  these  two  models. 


Table  7.  Missile  model  flight  conditions. 


Speed  (mph) 

Mach  no. 

q  (lbf/ff) 

60 

0.082 

9.77 

100 

0.136 

26.50 

Table  8.  Lambda  UCAV  flight  conditions. 


Speed  (mph) 

Mach  no. 

q  (lbf/ff  ) 

Re 

60 

0.082 

9.72 

3.5e5 

80 

0.109 

17.13 

4.6e5 

The  values  for  the  lift  coefficient,  Cl,  as  a  function  of  angle  of  attack,  a,  for  the 
bare  missile  are  shown  in  Figs.  22  and  23.  The  results  agree  somewhat  in  the  60  mph 
case  and  much  better  in  the  100  mph  case.  The  main  problem  is  that  even  at  its 
maximum  the  missile  is  only  producing  around  0.2  lbf  of  force,  depending  on  a.  This  is 
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too  small  a  number  to  expect  accurate  numbers  from  this  25-lb  balance,  which  was  also 
used  by  Dike.  One  possible  explanation  is  differences  in  how  the  model  was  installed. 

All  three  of  these  runs  were  accomplished  with  the  model  removed  in  between.  Different 
tare  runs  were  accomplished  for  each  run,  with  the  tare  possibly  having  small 
discrepancies.  Finally,  the  calibrations  done  for  this  study  and  the  Dike  [25]  study  were 
slightly  different.  All  of  these  errors  more  than  likely  combined  to  give  the  differences 
shown  in  the  curves.  The  missile  with  the  wing  installed  produced  significantly  more  lift 
than  the  bare  missile.  These  tests  also  prove  some  measure  of  repeatability  because  run 
one  was  accomplished  as  the  first  test  done,  and  run  two  was  accomplished  after  all  other 
runs  were  complete. 


alpha 


Figure  22.  Bare  missile  Cl  vs  alpha  -  60  mph. 
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Figure  23.  Bare  missile  Cl  vs  alpha  -  100  mph. 


In  order  to  try  and  get  a  better  measure  the  accuracy  of  the  balance  and  data 
reduction,  the  lambda  UCAV  model  was  tested  at  identical  conditions  to  the  ones  used  by 
In  [29].  The  results  can  be  seen  in  Figs.  24  and  25. 
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Figure  24.  My  lambda  UCAV  results. 


Figure  25.  In’s  lambda  UCAV  results  [In]. 
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My  results  match  up  very  well  to  the  Reed  results  shown  by  In  [29].  They  differ 
somewhat  in  the  values  of  the  other  runs,  but  the  shape  is  almost  identical.  The 
differences  are  probably  due  to  the  fact  that  In  used  a  100-lbf  balance,  where  a  25-lbf 
balance  was  used  in  this  experiment,  giving  slightly  different  values. 

Middle  Wing 

The  lift  coefficient  is  plotted  versus  angle  of  attack  in  Figs.  26  and  27.  In  these 
figures  the  sweep  angle,  X,  was  varied  from  0  to  90  deg  in  15  deg  increments.  For  all 
cases  S,  the  reference  area  used  to  calculate  Cl  was  the  planform  area,  and  was  kept 
constant  regardless  of  the  sweep  angle.  At  both  speeds  the  wing  stalled  at  X  =  0  and  15 
deg  at  around  a  of  4  deg,  but  then  recovered  causing  Cl  to  once  again  begin  to  rise.  This 
can  be  seen  in  the  large  drop  off  of  Cl.  The  wing  stalling  could  have  also  been  what 
caused  the  model  to  oscillate  up  and  down  in  the  original  100  mph  run  at  a  =  4  deg.  In 
this  configuration  the  X  =  45  deg  configuration  still  performs  well,  delaying  stall  until 
around  14  deg.  CL  maximum  occurs  at  X  =  0  deg. 

The  value  of  Cl  at  stall  at  X  =  0  was  slightly  higher  for  the  50  mph  case  than  for 
the  60  mph  case.  The  maximum  values  for  Cl  at  X  =  0  are  virtually  identical  for  50  and 
60  mph. 
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Figure  26.  Middle  wing  Cl  vs  alpha  -  50  mph. 
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Figure  27.  Middle  wing  Cl  vs  alpha  -  60  mph. 
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Cd  is  plotted  versus  a  in  Figs.  28  and  29  for  the  middle  wing  configuration. 

These  plots  look  like  typical  Cd  versus  a  plots.  The  drag  reaches  a  minimum  value  but 
never  goes  negative  and  then  steadily  increases  as  the  angle  of  attack  increases.  Also,  at 
high  angles  of  attack  the  drag  increases  as  X  decreases.  This  agrees  with  Tang  et  al.  [34] 
who  showed  that  for  an  oblique  wing  projectile  traveling  at  M  =  0.8,  the  drag  is  inversely 
proportional  to  sweep  angle.  One  possible  explanation  could  be  that  that  the  decrease  in 
friction  drag  caused  by  less  frontal  area  being  presented  to  the  air  stream  as  the  sweep 
angle  increases  is  larger  than  the  increase  in  pressure  drag  caused  by  the  fact  that  the 
wing  is  becoming  less  streamlined  as  it  is  swept,  although  no  pressure  measurements 
were  taken  in  this  study  to  verify  this.  One  puzzling  result  is  at  both  speeds  both  X  =  15 
and  30  deg  have  a  higher  Cd  than  X  =  0  deg  at  a  =  0  deg. 


-0  deg 
-15  deg 
30  deg 
45  deg 
-60  deg 
-75  deg 
-90  deg 


Figure  28.  Middle  wing  Cd  vs  alpha  -  50  mph. 
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Figure  29.  Middle  wing  Cd  vs  alpha  -  60  mph. 


Lastly,  Cm  is  plotted  as  a  function  of  a  for  the  middle  wing  configuration, 
measured  about  the  center  of  gravity.  This  can  be  seen  below  in  Figs.  30  and  3 1 . 
According  to  Nelson  [31],  in  order  for  an  aircraft  to  have  longitudinal  stability,  the  slope 
of  the  Cm  versus  alpha  curve  should  be  negative.  This  is  so  that  if  the  aircraft  encounters 
an  upward  gust  it  would  tend  to  want  to  respond  nose  down  therefore  returning  itself 
back  toward  equilibrium  [31].  With  the  exception  of  high  angles  of  attack  all  of  the 
slopes  are  negative  or  near  zero.  At  both  speeds  at  X  =  0  the  slope  is  very  negative 
indicting  very  good  pitch  stability.  Most  of  the  other  configurations  are  only  slightly 
negative  or  near  zero  so  the  tail  could  still  have  to  be  enlarged  or  moved  farther  back. 


37 


— 

0  -5  *  ^  """T 

alpha 


— ♦—  0  deg 
—1—15  deg 
30  deg 
45  deg 
-*-60  deg 
—•—75  deg 
— 1—90  deg 


Figure  30.  Middle  wing  Cm  vs  alpha  -  50  mph. 
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Figure  31.  Middle  wing  Cm  vs  alpha  -  60  mph. 


In-ground  effects 

The  missile  and  wing  were  run  upside  down  over  a  ground  plane  in  order  to 
simulate  release  from  an  aircraft.  This  puts  the  missile  in-ground  effect.  The  model  was 
run  at  two  different  speeds,  50  and  60  mph,  and  at  two  different  heights  with  h/b  ranging 
from  0.071  with  plane  2  at  X  =  0  to  1.94  with  plane  1  and  X  =  90.  The  h/b  is  different  for 
every  configuration  because  as  X  varies  so  does  b.  A  complete  list  of  all  h/b  values  can 
be  found  in  Table  6. 

The  values  of  Cl  versus  alpha  for  all  different  configurations  with  the  wing  in 
middle,  including  the  runs  without  the  ground  plane,  are  shown  below  in  Figs.  32  to  38. 
For  all  cases  plane  one  refers  to  the  ground  plane  farther  away  from  the  wing.  In  all 
cases  the  lift  goes  up  dramatically  as  the  ground  planes  get  closer  to  the  model.  This  was 
to  be  expected.  It  is  somewhat  hard  to  get  a  true  idea  of  what  is  occurring  because  the 
ground  plane  tests  had  to  be  limited  due  to  the  model  coming  in  contact  with  the  plane 
above  certain  angles  of  attack.  With  plane  one  installed  the  model  could  only  be  run  to  -5 
deg  angle  of  attack,  which  after  swapping  the  sign  becomes  +5  deg,  because  any  less 
angle  of  attack  would  have  caused  the  missile  nose  to  strike  the  ground  plane.  Similarly, 
with  plane  two  installed  the  angle  of  attack  could  not  be  lower  than  -2  deg.  It  should 
also  be  noted  that  for  X  =  90  deg  and  plane  one  installed  the  h/b  is  1 .94  which  is  outside 
of  the  normally  recognized  limit  for  IGE  of  h/b  =  1.  Although  this  is  the  case,  you  can 
still  see  a  clear  increase  in  lift  coefficient  with  ground  plane  one  installed  in  Fig.  38. 

For  all  cases  there  is  not  much  difference  between  the  data  at  50  and  60  mph.  At 
X  =  0  deg  the  difference  between  the  IGE  and  OGE  results  is  not  as  dramatic  as  it  was  for 
the  other  configurations,  although  it  did  stall  at  2  deg  IGE  as  opposed  to  4  deg  OGE.  The 
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smallest  increases  in  Cl  are  at  A,  =  0  and  90  deg,  with  the  largest  being  at  A  =  30,  45  and 
60  deg.  Also,  for  the  cases  where  there  is  a  large  difference  between  Cl  IGE  and  OGE, 
most  notably  A,  =  30,  45  and  60  deg,  the  increase  between  plane  one  and  plane  two  is 
nowhere  near  as  large  as  the  difference  between  OGE  and  plane  one.  The  IGE  curves  for 
0,  15,  30,  45,  and  60  deg  are  all  pretty  similar.  The  increase  in  Cl  is  then  caused  by  the 
differences  in  the  OGE  curves. 
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Figure  32.  Middle  wing  Cl  vs  alpha  -  k  =  0  deg  IGE. 
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Figure  33.  Middle  wing  CL  vs  alpha  -  k  =  15  deg  IGE. 
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Figure  34.  Middle  wing  Cl  vs  alpha  -  k  =  30  deg  IGE. 
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Figure  35.  Middle  wing  Cl  vs  alpha  -  k  =  45  deg  IGE. 
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Figure  36.  Middle  wing  Cl  vs  alpha  -  k  =  60  deg  IGE. 
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Figure  37.  Middle  wing  CL  vs  alpha  -  X,  =  75  deg  IGE. 


—♦—50 

mph 

—■—60 

mph 

—A— 50 

mph 

plane  1 

60 

mph 

plane  1 

— *—  50 

mph 

plane  2 

-0-60 

mph 

plane  2 

Figure  38.  Middle  wing  CL  vs  alpha  -  k  =  90  deg  IGE. 
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Next  the  drag  coefficient  variation  with  ground  plane  height  was  analyzed  for 
every  middle  wing  configuration.  The  plots  of  Cd  versus  a  are  shown  in  Figs.  39  to  45. 
These  plots  are  zoomed  in  on  the  area  where  the  ground  plane  measurements  were  taken 
to  get  a  better  picture  of  what  is  happening.  With  the  exception  of  X  =  75  and  90  deg,  the 
values  of  Cd  increase  as  the  ground  plane  gets  closer  as  long  as  a  is  positive.  This  was  a 
somewhat  expected  result  when  looking  at  the  Cl  curves  when  IGE.  A  major  contributor 
to  the  overall  drag  is  the  induced  drag,  also  known  as  drag  due  to  lift.  The  results  don’t 
match  very  well  with  the  results  predicted  using  McCormick’s  induced  drag  factor. 
According  to  his  formula,  any  h/b  value  of  less  than  around  0.5  should  give  a  lower 
induced  drag  [14].  In  this  study,  all  h/b  values  tested  gave  a  higher  induced  drag,  which 
does  agree  with  Lee  [18]  and  Curry  [20]  who  state  than  any  increase  in  drag  is  possible 
while  IGE. 
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Figure  39.  Middle  wing  CD  vs  alpha  -  k  =  0  deg  IGE. 
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Figure  40.  Middle  wing  Cd  vs  alpha  -  X  =  15  deg  IGE. 


45 


— ♦ —  50 

mph 

-■-60 

mph 

50 

mph 

plane  1 

60 

mph 

plane  1 

— *—  50 

mph 

plane  2 

o 

CD 

I] 

mph 

plane  2 

-10  -5  0  5  10  15  20 

alpha 


Figure  41.  Middle  wing  CD  vs  alpha  -  X,  =  30  deg  IGE. 
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Figure  42.  Middle  wing  CD  vs  alpha  -  k  =  45  deg  IGE. 
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Figure  43.  Middle  wing  Cd  vs  alpha  -  k  =  60  deg  IGE. 
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Figure  44.  Middle  wing  Cd  vs  alpha  -  k  =  75  deg  IGE. 
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Figure  45.  Middle  wing  CD  vs  alpha  -  A  =  90  deg  IGE. 

A  very  important  consideration  in  the  release  of  weapons  is  the  variation  of  Cm 
with  the  ground  plane  in  place.  This  will  dictate  if  the  missile  is  going  to  tend  to  want  to 
fly  back  up  into  the  aircraft  in  the  presence  of  a  small  disturbance.  The  Cm  versus  a 
curves  are  shown  in  Figs.  46  to  52.  The  model  shows  much  better  longitudinal  stability 
IGE  than  it  did  OGE.  The  only  configuration  of  concern  would  be  at  A,  =  0,  shown  in  Fig. 
46.  The  slope  of  the  line  is  virtually  zero.  Assuming  the  missile  would  be  released  at  A  = 
90  and  then  morph  into  a  A  =  0  configuration,  the  missile  should  be  well  clear  of  the 
aircraft  before  any  effects  would  be  seen  and  the  model  performed  very  well 
longitudinally  at  A  =  0  OGE.  All  of  the  other  configurations  had  a  nice,  negative  slope 
which  is  desirable  for  pitch  stability. 
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Figure  46.  Middle  wing  Cmvs  alpha  -  k  =  0  deg  IGE. 
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Figure  47.  Middle  wing  Cmvs  alpha  -  k  =  15  deg  IGE. 


49 


Cm  Cm 


— ♦ —  50  mph 
■  60  mph 

50  mph  plane  1 
60  mph  plane  1 
— *—  50  mph  plane  2 
•  60  mph  plane  2 


Figure  48.  Middle  wing  Cmvs  alpha  -  k  =  30  deg  IGE. 
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Figure  49.  Middle  wing  Cmvs  alpha  -  k  =  45  deg  IGE. 
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Figure  50.  Middle  wing  Cmvs  alpha  -  k  =  60  deg  IGE. 


Figure  51.  Middle  wing  Cmvs  alpha  -1  =  75  deg  IGE. 


Figure  52.  Middle  wing  Cmvs  alpha  -  X  =  90  deg  IGE. 
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Front  wing 

There  were  additional  attachment  points  at  the  front  of  the  missile  used  in  the 
joined  wing  experiments,  so  the  wing  was  mounted  at  that  point  and  the  tests  repeated. 
Figs.  53  and  54  show  Cl  as  a  function  of  a  for  50  and  60  mph.  Once  again  the  wing 
stalls,  but  this  time  at  an  angle  of  attack  of  2  deg,  and  this  configuration  does  not  recover 
as  quickly  as  when  the  wing  is  in  the  middle.  This  causes  a  lower  Cl  maximum  than  with 
the  middle  wing.  The  maximum  Cl  at  a  =  0  is  just  under  0.3  and  occurs  at  X  =  0  with  Cl 
decreasing  as  X  is  increased,  as  would  be  expected.  There  is  not  much  difference  between 
the  curves  for  50  and  60  mph.  At  both  speed  the  maximum  Cl  occurred  at  X  =  45  deg,  as 
opposed  to  occurring  at  X  =  0  deg  for  the  middle  wing  configuration.  Also,  the  maximum 
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value  of  Cl  is  only  about  one-half  of  the  maximum  value  for  the  middle  wing.  For  the 
sweep  angles  that  did  not  stall,  Cl  increased  as  the  sweep  angle  decreased. 
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Figure  53.  Front  wing  Cl  vs  alpha  -  50  mph. 
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Figure  54.  Front  wing  CL  vs  alpha  -  60  mph. 

The  drag  coefficient,  Cd,  is  plotted  as  a  function  of  a  in  Figs.  55  and  56.  All  the 
curves  have  the  same  basic  shape  with  the  exception  of  A.  =  0  and  15  deg  around  the  time 
when  stall  occurred  at  a  =  4  deg.  It  is  interesting  to  note  that  at  a  =  0,  Cd  values  are  in 
order  of  increasing  X  with  0  deg  having  the  lowest  value  and  90  deg  registering  the 
highest  drag.  These  values  switch  at  higher  angles  of  attack  with  X  =  0  having  the  highest 
value  at  a  =  15  and  X  =  90  having  the  lowest  value  of  Cd-  You  can  also  see  that  at  a  =  0 
deg  the  drag  on  the  X  =  0,  15,  30  and  45  deg  configurations  is  essentially  zero.  The  Cd 
for  X  =  0,  30,  45,  and  60  deg  drops  below  0,  which  is  also  puzzling.  This  could  partially 
be  due  to  the  disturbed  airflow  due  to  the  nose  causing  a  suction  effect.  Another  anomaly 
caused  by  the  interaction  between  the  wing  on  the  front  of  the  missile  and  the  nose  could 
also  be  the  cause.  Pressure  measurements  and  possibly  flow  visualization  are  needed  to 
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accurately  diagnose  the  cause  of  these  drag  curves.  The  shape  of  the  curves  also  differs 
from  parabolic  shape  normally  seen  in  Cd  versus  alpha  plots. 
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Figure  55.  Front  wing  CD  vs  alpha  -  50  mph. 
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Figure  56.  Front  wing  Cd  vs  alpha  -  60  mph. 
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The  next  parameter  to  be  examined  is  the  pitch  moment  coefficient,  Cm, 
measured  about  the  center  of  gravity,  CG,  of  the  model.  Cm  as  a  function  of  alpha  is 
shown  in  Figs.  57  and  58.  This  configuration  provides  much  worse  pitch  stability  than 
the  middle  wing  configuration.  This  is  primarily  due  to  the  placement  of  the  wing.  With 
the  wing  located  near  the  CG,  the  moment  arm  is  greatly  reduced  producing  a  much 
smaller  pitch  up  moment  due  to  the  lift  created  by  the  wing.  This  smaller  moment  can  be 
much  more  easily  counteracted  by  the  small  tail.  The  only  configuration  that  is  truly 
longitudinally  stable  is  at  X  =  90  deg  and  even  this  configuration  is  only  very  slightly 
stable.  The  rest  of  the  configurations  begin  longitudinally  unstable,  although  the  X  =  0 
and  15  deg  configurations  do  become  slightly  stable  after  stall.  There  is  also  very  little 
difference  between  the  50  and  60  mph  data. 

This  missile  was  designed  without  the  wing  so  therefore  has  a  relatively  small  tail 
when  compared  to  the  wing.  The  tail  is  too  small  to  overcome  the  upward  pitching 
moment  created  by  the  wing  producing  lift.  If  this  model  were  to  fly  the  tail  would  have 
to  be  redesigned  and  made  much  bigger  to  achieve  longitudinal  stability. 
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Figure  57.  Front  wing  Cm  vs  alpha  -  50  mph. 
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Figure  58.  Front  wing  Cm  vs  alpha  -  60  mph. 
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The  IGE  plots  for  the  front  wing  configuration  can  be  seen  in  Appendix  B. 
Overall  the  increase  in  Cl  IGE  for  the  front  wing,  shown  in  Figs.  59  to  65,  was  not  as 
large  as  for  the  middle  wing  and  in  the  case  of  X  =  60,  75,  and  90  deg  it  actually  went 
down.  The  Co  increase  IGE,  shown  in  Figs.  66  to  72,  was  also  much  higher  than  the 
increase  for  the  middle  wing  for  most  configurations.  The  IGE  Co  values  were  similar 
for  the  front  and  middle  wings,  but  the  OGE  values  for  the  front  wing  were  so  low  that 
the  net  increase  from  OGE  to  IGE  was  higher.  Also,  with  the  ground  planes  installed  the 
drag  never  goes  negative  like  it  does  OGE.  The  pitch  moment  performance  is  bad  IGE, 
shown  in  Figs.  73  to  79,  with  the  front  wing,  like  it  was  OGE. 

Uncertainty  Analysis 

The  uncertainty  of  a  function  of  several  variables  can  be  found  using  the  partial 
derivatives.  For  instance,  if  y  =  f(xi,X2,. . .)  the  uncertainty  dy  can  be  found  using 
equation  (3)  [33]. 

dy  =  )2  +  #  A,  )2  +  ...  +  #&,  )2  (3) 

y  ax  j  ax  2  ax  n 

For  this  analysis,  the  uncertainties  of  the  individual  measurements  were  given  in 
manufacturer  specifications  for  the  devices  used  to  make  the  measurements  and  are 
shown  in  Table  9.  The  wing  area  and  angle  of  attack  are  assumed  to  be  accurate. 
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Table  9.  Measurement  uncertainties. 


Measurement 

Uncertainty 

N 

0.0625  lbf 

A 

0.075  lbf 

V 

0.25  mph 

P 

0.01  psi 

The  uncertainties  were  then  calculated  for  the  50  mph,  middle  wing  configuration 
for  the  lift  and  drag  coefficients.  The  maximum  uncertainty  for  Cl  was  0.038  and  for  Cd 
was  0.045.  The  maximum  value  of  Cl  for  this  test  was  0.75  giving  an  uncertainty  of 
approximately  5%.  For  this  test  the  maximum  value  of  Co  was  0.41  giving  an 
uncertainty  of  1 1%. 
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Conclusions  and  Recommendations 


Section  1  -  Conclusions 

The  model  was  tested  OGE  with  the  wing  installed  both  on  the  front  of  the  missile 
and  in  the  middle  and  the  sweep  angle  varying  from  0  to  90  deg.  The  wing  stalled  at  an 
angle  of  attack  of  either  2  or  4  deg  depending  on  the  wing  placement.  The  middle  wing 
configuration  produced  more  lift,  but  also  significantly  more  drag  than  the  front  wing. 

The  validity  of  the  drag  calculations  is  up  for  debate  due  to  the  fact  that  in  certain 
configurations,  a  negative  drag  was  measured.  In  order  to  get  a  true  understanding  of 
what  is  happening,  additional  testing  will  need  to  be  accomplished  using  pressure 
transducers  and  some  sort  of  flow  visualization.  The  model  also  showed  good 
longitudinal  stability  in  the  middle  wing  configuration,  but  poor  stability  with  the  front 
wing.  This  is  caused  by  the  tail  not  being  large  enough  to  overcome  the  increased 
moment  arm  between  the  lift  and  the  CG. 

The  model  was  also  tested  inverted  over  the  ground  plane  to  simulate  dropping 
after  release  from  an  aircraft.  With  the  wing  installed  in  the  middle  the  model  exhibited  a 
higher  Cl  but  also  a  higher  Co  than  in  identical  conditions  OGE.  The  higher  Cl  was 
expected,  but  the  increase  in  CD  was  somewhat  unexpected  even  though  several  studies 
have  shown  that  certain  aircraft  do  experience  higher  drag  1GE  than  OGE.  The  model 
also  performed  well  longitudinally  1GE  at  all  sweep  angles  except  for  0  deg.  This 
wouldn’t  cause  a  problem  because  assuming  the  missile  is  dropped  from  the  aircraft  at  X 
=  90  deg  and  then  swept  to  the  optimal  angle,  if  it  did  sweep  to  X  =  0  deg,  the  missile 
would  be  well  clear  of  the  aircraft  where  it  is  stable  OGE. 
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The  IGE  tests  with  the  wing  in  the  front  showed  an  increase  in  lift  for  most 
configurations.  There  were  a  couple  of  instances  where  the  lift  was  less  IGE  than  OGE. 
The  drag  increase  for  this  configuration  was  significant,  with  the  drag  curves  looking 
more  like  traditional  Cd  versus  alpha  plots.  The  drag  also  never  goes  to  zero  IGE.  Much 
like  the  OGE,  the  front  wing  setup  performed  poorly  IGE  with  regards  to  longitudinal 
stability. 

It  was  tough  to  get  a  true  idea  of  the  full  ground  effects  on  this  model  because  of 
the  small  range  of  angles  of  attack  that  could  be  tested  due  to  model-ground  plane 
interferences.  It  would  also  be  helpful  to  run  the  missile  inverted  without  the  ground 
plane  to  verify  the  approach  used  in  analyzing  the  data,  and  also  to  rule  out  any  other 
effects.  Since  the  true  effects  of  oblique  wings  are  seen  at  transonic  and  supersonic 
speeds,  this  study  needs  to  be  repeated  at  those  speeds  because  these  results  would  not 
apply  to  the  more  complex  aerodynamics  in  these  regimes. 

Section  2  -  Recommendations 

This  study  only  provides  a  first  look  at  ground  effects  for  a  missile  with  an  oblique  wing 
installed.  Based  on  the  results  and  findings  of  this  study,  the  following  are 
recommendations  for  further  research: 

-  redesign  the  airfoil  to  avoid  the  wing  stall 

-  strengthen  the  wing  along  its  span  in  order  to  be  able  to  test  at  higher  Re 

-  add  some  method  of  flow  visualization  to  see  the  span  wise  airflow 

-  accomplish  a  CFD  analysis  to  validate  results 

-  re-accomplish  these  tests  at  both  transonic  and  supersonic  speeds 
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run  front  wing  tests  with  pressure  transducers  and  flow  visualization  to  try  and  diagnose 
cause  of  negative  drag 

design  a  smaller  model  to  allow  for  more  ground  plane  heights  and  angles  of  attack  to 
be  tested 
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Appendix  A:  Sample  Calculations 


Data  Reduction 

Test  condition:  Missile  with  no  wing  installed 

U  =  105.248  mph  =  154.364  ft/s 
a  =  12.348  deg 

Test  room  conditions: 

T  =  534.07  R 
P  =  14.278  psi 
R  =  1716  ft-lb/slug-R 
p  =  .372x1  O'6  slug/ft-s 
y  =  1.4 


S  =  ref  area  =  0.4  ft“ 

p  lb 

p  = - =  .00224 — y 

RT  ft 3 

q  =  -pU2  =26.68  l\ 
2  ft 2 


=  -N//RT  =1132.72  — 
sec 


Corrected  forces  obtained  from  MATLAB 
A  =.1024  lbf 
Y  =  -.004  lbf 
N  =  -1.2158  lbf 
1  =  -.0035  in-lbf 
m  =  1.9175  in-lbf 
n  =  0.008  in-lbf 

Angles  used 

0  =  a  =  12.348  deg 
\\i  =  0  deg 

Now  to  convert  from  body  axes  to  wind  axes 


~D~ 

A  cos  9  cos  y/  +  Y  sin  y/  +  N  sin  6  cos  y/ 

S 

= 

-A  sin  y/  cos  0  +  Y  cos  y/  -  N  sin  9  sin  yr 

L 

-A  sin  9  +  N  cos  9 
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D  =  .3062  lbf 
S  =  .004  lbf 
L  =  .6636  lbf 

Add  blockage  corrections  [32] 

^„  =  K^Body: =.005489 
c 

Ucorr  =  U  (1  +  s)  =  155.369  ft/s  =  105.934  mph 
cjcorr  =  q  (1  +  s)  =  26.974  lbf/ ft2 


C,  =  — ^ — =.061 

1 

^Vcon-S 

CD  =  - — =  .028 

^  Qcorr^ 


M 


U, 


=  .137 


a 


Re  =  pU'0”1 
A 


=  2.22xl06 


Error  Analysis 


CL  = 


L 

-pV2S 

2 


and 

L  =  N  cos  a  -  A  sin  a 


so, 
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N  cos  a  -  A  sin  a 

-pV2S 

2 


The  uncertainty  is  given  by: 


With  the  individual  partial  derivatives: 


ac, 

ON 

ac, 

8A 


cos  a 

-pV2S 

2 

-sin  a 


0.591 

-0.108 


ac,  -l(iVcosa-^sina)  „ 

— -  =  — - — - -=  -0.014 

8P  ip’V’S 

2 

dCL  _  -2(Ncosa  -Asm  a)  _ 

W=  i 

2 

Multiplying  by  the  uncertainties  for  each  measurement  given  in  Table  9  gives  dCL  = 
0.038. 
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Appendix  B.  Additional  Plots 
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Figure  59.  Front  wing  ClVS  alpha  -  k  =  0  deg  IGE. 
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Figure  60.  Front  wing  CLvs  alpha  -  k  =  15  deg  IGE. 
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Figure  61.  Front  wing  ClVS  alpha  -  k  =  30  deg  IGE. 
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Figure  62.  Front  wing  ClVS  alpha  -  k  =  45  deg  IGE. 
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Figure  63.  Front  wing  CLvs  alpha  -  k  =  60  deg  IGE. 
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Figure  64.  Front  wing  CLvs  alpha  -1  =  75  deg  IGE. 
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Figure  65.  Front  wing  CLvs  alpha  -  k  =  90  deg  IGE. 
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Figure  66.  Front  wing  CdVS  alpha  -  k  =  0  deg  IGE. 
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Figure  67.  Front  wing  CDvs  alpha  - 1  =  15  deg  IGE. 
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Figure  68.  Front  wing  CdVS  alpha  -  k  =  30  deg  IGE. 
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Figure  69.  Front  wing  Covs  alpha  -  k  =  45  deg  IGE. 
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Figure  70.  Front  wing  CdVS  alpha  -  k  =  60  deg  IGE. 
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Figure  71.  Front  wing  Covs  alpha  -k  =  75  deg  IGE. 


— ♦ — 50 

mph 

■  60 

mph 

. 

o 

mph  plane  1 

60 

mph  plane  1 

i 

o 

mph  plane  2 

•  60 

mph  plane  2 

-10  -5  0  5  10  15  20 

alpha 


Figure  72.  Front  wing  Covs  alpha  -  k  =  90  deg  IGE. 
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Figure  73.  Front  wing  Cmvs  alpha  -  k  =  0  deg  IGE. 
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Figure  74.  Front  wing  Cmvs  alpha  -  k  =  15  deg  IGE. 
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Figure  75.  Front  wing  Cmvs  alpha  - 1  =  30  deg  IGE. 
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Figure  76.  Front  wing  Cmvs  alpha  -  k  =  45  deg  IGE. 
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Figure  77.  Front  wing  Cmvs  alpha  -  k  =  60  deg  IGE. 
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Figure  78.  Front  wing  Cmvs  alpha  -X  =  75  deg  IGE. 
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Figure  79.  Front  wing  Cmvs  alpha  -  k  =  90  deg  IGE. 
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Appendix  C:  MATLAB  Data  Reduction  Program 


l  k  ~k  -k  -k 
l  -k  ~k  ~k  -k 
t-k-k-k-k 
t-k-k-k-k 
L  -k  ~k  -k  ~k 
'^-k  ~k  -k  -k 
l  ~k  ~k  ★ 
t-k-k-k 
L  Vr  ~k  ~k  ~k 
L  ~k  ~k  ~k  ~k 
l-k-kicic 
l-k-k-k-k 


'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

******  Lt .  Dave  Gebbie  &  Capt  Anthony  DeLuca  ******************* 
******  Adapted  fpr  the  Balance  AFIT  1  by  Lt.  Rivera  Parga********* 

******  Calculation  of  Lift,  Drag,  Moments  ******************* 

******  Re-adapted  by  Troy  Leveron,  ENS,  USNR  ******************* 
Re-adapted  by  Brett  Jones,  ENS,  USNR  for  UCAV  Ground  Effects  Test**** 
Re-adapted  by  Won  In,  Capt,  USAF  for  UCAV  Ground  Effects  Test******** 
Re-adapted  by  Matt  Dillsaver,  lLt,  USAF  for  Oblique  Wing  Test******** 
************** *cal CUla tion  of  Lift,  Drag,  Moments********************* 

•k'k-k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k-k'k'k'k'k'k'k'k'k'k'k 

^^^^^^^^^^iciir^^iiciic^^i^'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k-k'k'k'k'k'k'k 


%This  Code  will  transfer  measured  Forces  and  Moments  on  the  balance  to  Wind 
% (earth) centered  frame  of  reference  by  correctiing  for  tare  eff ects , balance 
%interactions ,  and  wind  tunnel  irregularities,  then  plot  lift  and  drag 
%coef ficients  in  as  functions  of  AoA 

clear;  clc;  format  long; 

^^kcicic^^i^iir^kcici^^^ici^^^icic^^i^ic^'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

%INPUT  DECK 

%FIRST  FILL  THE  FOLLOWING  INFORMATION 

9 -i^kri^ic^^i^i^^^ic'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 


Masskg  =  3.61;  %KGs 

Troom  =  73.1  +  459.67;  %deg  R 
P~barro  =  29.2262  *  0.4911541; 


****Changed  for  each  day  of  testing**** 

%Psi  *Changed  for  each  day  of  testing**** 


%  INPUT  DATA  FILE  AND  INPUT  DATA  TARE  FILE 
load  MATT2OW0tareA_4tol5B0.txt; 
TareFile=MATT2OW0tareA_4tol 5B0 ( : , 1: 9) ; 
load  MATT20W0M50A_4tol5B0.txt; 
DataFile=MATT2OW0M50A  4tol5B0  ( : , 1: 9) ; 


%  Raw  tare  file  to  be  read  in. 
%  Raw  data  file  to  be  read  in: 


%Offset  distances  from  the  Mounting  Block  to  the  Model  C.G.  (inches) 


Y_cmb 
X_cmb 
Z  cmb 


0.0; 

1.5625; 

0.0; 


%inches  (from  origin  @  balance  center  w /  +  right) 
%inches  (from  origin  @  balance  center  w/  +  down) 


%Required  for  the  Solid  Body  blockage  corrections  due  to  wing  and  fuselage 
Body_Volume  =  .4;  %ftA3: 

Wing  Area  =  35/12A2;  %ftA2 

%I . -  Room  Conditions  and  Model  Specifics 
%  Units  are  in  ft,  sec,  lbm,  psf,  rankine,  fps 


Mass=  (Masskg  *  1000)  *  .0022046;  %  lbm 

Gas_Const  =  1716;  %  f t-lbf /slug-R 

Density  =  (P_barro  *  144)/ (1716  *  T  room) ;  %lbm/ftA3  or  lbf-sA2/ftA4 
Root_Chord  =  2/12  %ft 
Span  =  17.5/12  %ft 
Aspect_Ratio  =  SpanA2/Wing_Area; 
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Kinematic_Viscosity  =  .372e-6;  %slug/ft-s 
Speed_of_Sound  =  sqrt(1.4  *  Troom  *  Gas_Const) ;  %fps 


oDistances  between  sensors  (inches)  to  calculate  moments 


D1 


(2.10  /  2);  D2  =  Dl;  D3  =  (1.7  /  2);  D4  =  D3 ; 


%inches 


^'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

%II. -Solid  Body  Blockage  corrections  due  to  wing  and  fuselage  (Pope  pg  369) 

^■k'k'k'k'k-kic-k-k'k-kic'k-k'kic'k'k'k'k'k-k-k-k-k'k-k-k'k-k-k-k-k'k-k-k'k-k-k-k-k'k-k-k'k-k-k-k-k'k-k-k'k'k-k-k-k'k-k-k-k-k-k-k-k-k-k-k'k'k-k-k-k-k 

K_1  =  1.12;  %  From  fig.  10.2  assuming  4  digit  airfoil  and  t/c=.15 

delta  =  0.122;  %  boundary  correction  factor  (b/B)  Ch  10 

Tau_l  =  0.89;  %  factor  from  pg  369  fun.  of  tunnel  shape  and  b/B 

X_Section  =  (32/12) *(  41/12)  ;  %ftA2 

Wing  Volume  =  29.6861/12A3  %  ftA3 

K_3  =  0.89; 

Epsilon_sb_w  =  (K_l*Tau_l*Wing_Volume)  /  X_SectionA (3/2 )  ; 
Epsilon_tunnel_correction=l ; 

Epsilon_sb_b  =  (K_3*Tau_l*Body_Volume)  /  X_SectionA (3/2)  ; 

Epsilon_tot  =  Epsilon_sb_w  +  Epsilon_sb_b+ (Epsilon_tunnel_correction-l ) ; 


Q-'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

%VI.  -CORRECT  FORCES  AND  MOMENTS  FOR  BALANCE  INTERACTIONS  (BODY  AXIS) 

%Balance  Interactions  with  off  axis  elements  for  the  25  lb  balance 

%Using  average  of  the  25  lb  calibration  runs  for  N1  &  N2  and  the 

%10  lb  calibration  for  SI,  S2  &  A  then  normalizing  by  the  actual 

%sensor  (Nl,  N2,...)  in  question.  The  sensor  sequence  in  each  row  vector  is: 

% [N1  N2  SI  S2  A1  L2 ] 


N1_I  =  ([1  -.010272  .01244  -.000593  .13758  .014914] ) ; 

Nil  =  N1_I (1,1); 

N2_I  =  ([-.023  1  .00079  -.000296  .019354  .14516]); 

N22  =  N2_I (1,2); 

S1_I  =  ([.002864  .00593  1  .002963  .089481  .000099]);  Sll  =  S1_I(1,3); 

S2_I  =  ([-.001576  .003939  -.026194  1  .028065  .007287]); 

S22  =  S2_I (1,4); 

A_I  =  ([.209628  -.21866  .008038  -.007135  1  -.002077]); 

All  =  A_I (1,5) ; 

L_I  =  ([-.253931  -.239861  -.007284  .019533  .025824  1] ) ; 

LI 1  =  LI (1,6) ; 


Nl 

normalized  = 

(Nl  I) 

.*  [Nil  N22 

Sll  S22  All  LI 1 ] 

•  A  (-D  ; 

N2 

normalized  = 

(N2  I) 

.*  [Nil  N22 

Sll  S22  All  LI 1 ] . 

A  (-D  ; 

SI 

normalized  = 

(SI  I) 

. *  [Nil  N22 

Sll  S22  All  LI 1 ] 

■  A  (-D  ; 

S2 

normalized  = 

( S 2  I) 

. *  [Nil  N22 

Sll  S22  All  LI 1 ] 

•  A  (-D  ; 

A 

normalized  = 

(A  I) 

.*  [Nil  N22 

Sll  S22  All  LI 1 ] . 

A  (-D  ; 

L 

normalized  = 

(L  I) 

.*  [Nil  N22 

Sll  S22  All  LI 1 ] . 

A  (-D  ; 

Interactions  Kij  =  [Nl 

normalized ' 

N2  normalized'  SI 

normalized ' 

S2  normalized' 


A_normalized'  L_normalized' ] ; 


%Load  the  static  tare  data  for  the  alpha  sweep  w/o  the  wind  ,  separate  each 
%force  from  the  file,  and  fit  a  4th  order  poly  as  an  x-y  plot  (AoA  vs. 
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%Force)  for  each  of  the  6  force  sensors. 

SJj' k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 


%load  tarel.txt;  %raw  tare  data  file  to  be  read  in 
FILE=TareFile  ( : , 1 : 9) ; 


j=i; 

k=l; 

L=length (FILE) ; 

for  i=l:L 
if  i~=L 

next 

NEXT=i+l; 

next  row 

VALUE2=FILE (NEXT, 1) ; 
else  if  i==L 
VALUE 2 =50; 

end 

end 

A ( j ,  : ) =FILE  (i,  :  )  ; 

VALUE1=FILE  (i, 1) ; 

FILE 

if  VALUE 1 ==VALUE 2 

row 

j=j+i; 

else  if  VALUE 1 ~ =VALUE 2 
check 

if  length (A ( : , 1 ) ) <5 
due  to  angle  change 

j=i; 

clear  A; 

else  if  length (A ( : , 1 ) ) >5 
C=length (A ( : , 1) ) ; 
for  m=l : 9 


end 

j=i; 

k=k+l ; 
clear  A 

end 

end 

end 

end 

end 


%Run  for  all  data  points  #  of  rows 
%if  current  row  is  not  last  row,  go  to 

%set  next  equal  to  the  value  of  the 

%set  value2  as  next  row  column  1 
%unless  the  it  is  the  last  value 
%value2  set  to  50  to  end  the  sequence 


%set  row  j  of  A  equal  to  row  i  of  FILE 
%set  valuel  equal  to  row  i  column  1  of 

%if  valuel  equals  value2,  go  to  next 


%if  valuel  and  value2  are  different 

%if  less  than  20  values,  ignored 


%if  more  than  20  values 
%find  length  of  A 
%Average  all  rows  of  the  like 

(4:C,m));  %disregarding  first  10  for 


if  B (k-1, 1) <B ( ( k-2 ) , 1) 
B=B  ( 1  :  (k-2) ,  :) 

end 


tare= [B] 

%******************************End  of  Inserted  Code************************ 
[row, col]  =  size (tare); 


for  k  =  1 : row 

tare(k,l).*  (pi/180); 
tare (k, 4) ; 
tare ( k, 5)  ; 
tare (k, 7) ; 


theta_tare  (k,  : ,  :  ) 
Nl^tare (k, : ,  :  ) 
N2_tare (k,  :  ,  :  ) 
Sl_tare (k, : , : ) 
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=  tare (k, 8)  ; 
=  tare ( k, 6)  ; 
=  tare (k, 9)  ; 


S2_tare  (k,  :  ,  :  ) 

A_tare (k,  : ,  :  ) 

L_tare ( k, : , : ) 

end 

Nl_poly  =  polyfit (theta^tare, Nl_tare, 4) ; 

N2_poly  =  polyfit (theta_tare, N2_tare, 4 ) ; 

Sl_poly  =  polyfit (theta_tare, Sl_tare, 4 ) ; 

S2_poly  =  polyfit (theta^tare, S2_tare, 4) ; 

A_poly  =  polyfit (theta^tare, A_tare,  4 )  ; 

L_poly  =  polyfit (theta_tare, L_tare, 4 )  ; 

clear  ('B', 'C' , 'D' , 'L') 

%Load  the  specific  test  run  files,  subtract  the  effect  of  the  static 
%weight  with  the  tare  polynominals  above,  and  correct  for  cross  diagonal 
%sensor  interactions. 

FILE=DataFile (:,:);  %Same  as  Above 


j=i; 

k=l; 

L=length (FILE) ; 


for  i=l:L 
if  i~=L 

next 

NEXT=i+l; 

next  row 

VALUE2=FILE (NEXT, 1) ; 
else  if  i==L 
VALUE 2 =5 0; 

end 

end 

A ( j ,  :)=FILE  (i,  :  )  ; 

VALUE 1=FILE  (i, 1) ; 

FILE 

if  VALUE 1 ==VALUE 2 


%Run  for  all  data  points  #  of  rows 
%if  current  row  is  not  last  row,  go  to 

%set  next  equal  to  the  value  of  the 

%set  value2  as  next  row  column  1 
%unless  the  it  is  the  last  value 
%value2  set  to  50  to  end  the  sequence 


%set  row  j  of  A  equal  to  row  i  of  FILE 
%set  valuel  equal  to  row  i  column  1  of 

%if  valuel  equals  value2,  go  to  next 


row 

j=j+i; 

else  if  VALUE 1 ~  =VALUE 2 


%if  valuel  and  value2  are  different 


check 

if  length (A ( : , 1 ) ) <5 
due  to  angle  change 

j=i; 

clear  A; 

else  if  length (A ( : , 1 ) ) >5 
C=length (A ( : , 1) ) ; 
for  m=l : 9 

values  in  A 


vibrations 


B  ( k, m)  =mean  (A  ( 4  :  C, m)  )  ; 


%if  less  than  20  values,  ignored 


%if  more  than  20  values 
%find  length  of  A 
%Average  all  rows  of  the  like 

%disregarding  first  10  for 


end 


j=i; 
k=k+l; 
clear  A 


end 

end 


end 
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end 


end 

if  B (k-1,  1) <B ( ( k-2 ) , 1) 
B=B  ( 1  :  (k-2) ,  :) 

end 


sample_data= [B] 

%******************************End  of  Inserted 
Code*  *  *************************** 


[row2,col2]  =  size ( sample_data) ; 
for  i  =  1 : row2 


%Angles  of  the  model  during  test  runs  (Roll,  Pitch  {AoA},  Yaw  {Beta}) : 


phi  =  0 ; 

theta (i,:)  =  sample_data (i , 1 )  .*  (pi/180); 

si(i, :)  =  sample_data (i , 2 )  .*  (pi/180); 

Wind_Speed (i, ; )  =  sample_data (i , 3 )  .*  (5280/3600); 


%Flight  Parameters  (Re#,  Ma#,  Dynamic  Pressure) : 

q  =  (.5  *  Density)  .  *  Wind_Speed. A2 ; 
q_Corrected  =  q  .*  (1  +  Epsilon_tot ) A2 ; 

Wind_Speed_Corrected  =  Wind_Speed  .  *  (1  +  Epsilon_tot) ; 
Wind_Speed_Corrected_jnph  =  Wind_Speed_Corrected  .*  (3600/5280); 

Mach_Number  =  Wind_Speed_Corrected  ./  Speed_of_Sound;  %NonDimensional 

Reynolds_Number  =  ((Density  *  Root_Chord)  . *  Wind_Speed_Corrected)  ./ 
Kinematic_Viscosity ;  %NonDimensional 

Flight_Parameters  =  [Mach^Number  Reynolds_Number  q_Corrected] ; 

%individual  forces  for  each  sensor: 


%radians 

%radians 

%fps 


%lbf/ftA2 
%lbf /ft A2 
%fps 


INEW  NOTATION 


Nl_test (i, 
N2_test (i , 
Sl_test (i, 
S2_test (i, 
A_test (i,  : ,  : ) 
L  test  (i,  :  ,  :  ) 


)  =  sample_data (i , 4 ) ; 
)  =  sample_data (i , 5 ) ; 
)  =  sample^data (i , 7 ) ; 
)  =  sample_data (i , 8 ) ; 

=  sample^data (i, 6) ; 
=  sample_data (i , 9 ) ; 


k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

iv .  Subtract  the  effect  of  static  weight  w/  tare  polynomials  above 


%Evaluating  the  actual  test  theta  angle  (AoA)  in  the  tare  polynominal  to 
%determine  the  tare  values  for  the  angles  tested  in  each  run. 


Nl_eval  =  polyval (Nl_poly, theta) ; 
N2_eval  =  polyval (N2_poly, theta) ; 
Sl_eval  =  polyval (Sl_poly, theta) ; 
S2_eval  =  polyval (S2_poly, theta) ; 
A_eval  =  polyval (A_poly, theta) ; 
L_eval  =  polyval (L  poly, theta); 


%The  Time-Averaged  (raw)  forces  Nl,  N2,  SI,  S2,  A,  L  measurd  in  the  wind 
%tunnel  (body  axis)  with  the  tare  effect  of  the  weight  subtracted  off. 


Nl_resolved  =  Nl  test  -  (Nl_eval) ; 
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N2_resolved  =  N2_test  -  (N2_eval) ; 

Sl_resolved  =  Sl_test  -  (Sl_eval) ; 

S2_resolved  =  S2_test  -  (S2_eval) ; 

A_resolved  =  A_test  -  (A_eval); 

L_resolved  =  L_test  -  (L_eval); 

Forces_minus_tare  =  [Nl_resolved  N2_resolved  Sl_resolved  S2_resolved  A_resolved 
L_resolved] ' ; 

%Forces  Nl,  N2,  SI,  S2,  A,  &  L  corrected  for  the  balance  interactions  (body 
axis) 

Corrected_Data  =  ( inv ( Interactions_Ki j )  *  Forces_minus_tare)  ; 

S^'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

%Calculation  of  the  Axial,  Side,  &  Normal  Forces  from  the  corrected  balance 
%forces  in  the  Body  Axis  reference  frame 

^  ■ k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

Forces_b ( : , i)  =  [Corrected_Data ( 5 , i ) ;  Corrected_Data (3, i)  + 

Corrected_Data (4, i) ;  Corrected_Data ( 1 , i )  +  Corrected_Data (2 , i ) ] ; 

%Calculation  of  the  Drag,  Side,  &  Lift  Forces  in  the  Wind  Axis  reference 
% frame 


Forces_w  = 

[Forces_b (1 , : ) . *cos (theta ' ) . *cos ( si ' ) +Forces_b (2, :) . *sin (si ' ) +Forces_b (3, : ) . *si 
n  (theta ' )  . *cos  (si ' ) ; 

-Forces_b (1,  :)  . *sin (si ' )  . *cos (theta ' ) +Forces_b (2 , : )  . *cos  (si ' ) - 
Forces_b (3,  : )  . *sin  (theta ' )  . *sin  (si ' ) ; 

-Forces^b ( 1 ,  : )  . *sin (theta ' ) +Forces_b (3 ,  : )  . *cos  (theta ' ) ] ; 

%First  entry  is  the  moments  calculated  by  the  balance  or  direct  calculation 
%in  the  Body  Reference  Frame.  Balance  measures  Roll  (1),  Yaw  is  about  the 
%z-axis  (n) ,  and  Pitch  is  about  the  y-axis  (m) .  Distances  from  strain 
%gages  to  C.G.  are  in  INCHES.  Moments  are  in-lbf 

m  =  Corrected_Data (1, i)  *  D1  -  Corrected^Data (2 , i )  *  D2 ; 
n  =  Corrected_Data (3, i)  *  D3  -  Corrected_Data ( 4 , i )  *  D4; 

Moments_b ( : , i )  =  [Corrected_Data ( 6, i) ;  m;  n] ; 

%Second  entry  is  the  conversion  from  the  "Balance  Centeric"  moments  to  the 
%Wind  Reference  monments  with  respect  to  the  Balance  Center  (be) 

Moments_w_bc  =  [Moments_b ( 1 ,:). *cos (theta '). *cos ( si ') - 
Moment s_b (2,  :)  . *sin (si ' ) +Moments_b (3,  : )  . *sin (theta ' )  . *cos  (si ' ) ; 

Moment s_b (1,  :)  . *sin  (si ' )  . *cos (theta ' ) +Moments_b (2,:).*cos(si') +Moments_b (3,  : )  . * 
sin  (theta ' )  . *sin  (si ' ) ; 

-Moments_b (1 ,  : )  . *sin (theta ' ) +Moments_b  (3, : )  . *cos  (theta ' ) ] ; 

%Finally,  the  balance  centered  moments  are  converted  to  moments  about  the 
%Model ' s  Center  of  Mass  (cm)  or  Center  of  Gravity  (CG) 


cgdist=sqrt ( (X_cmb) A2+ ( Z_cmb) A2 ) ;  %Obtaining  the  direct  distance  between  the 

%center  of  the  balance  and  the  center  of 


mass 

w=atan ( -Z_cmb/X_cmb ) ;  %Obtaining  the  angle  between  cgdist  and  the  x  axes  at 
zero  angle  of  attack 


X_cm(i,  : ) =  cos (theta  (i,  : ) +w) *cos (si (i, : ) ) * (cgdist) ; 

Y_cm(i,:)  =  Y_cmb  +  X_cm (i,  : ) *tan  (si (i,  : ) ) ;  %  appropriate  for  very 

small  y_cmb  and  reasonable  si 
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Z  cm(i. 

:)=  -sin  (theta (i, :) +w) * (cgdist) 

i  ; 

Moments 

w  eg  u  =  [Moments  w  bc(l, :] 

)  + 

Z  cm ( i , 

Y  cm ( i , 

: )  ; 

Z  cm(i. 

Moments  w  be (2,;) 

:  )  ; 

)  - 

X  cm ( i , 

X  cm ( i , 

Moments  w  be (3,  :  [ 

:  )  ] ; 

)  - 

Y  cm ( i , 

*Forces_w (2, : ) 
*Forces_w (3, : ) 
*Forces_w ( 1 , : ) 


+  Forces^w (3, : ) * 
+  Forces^w ( 1 , :  )  * 
-  Forces  w (2, : ) * 


SJj" k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

%Calculation  of  the  actual  Lift  and  Drag  nondimensional  Coefficients, 
%uncorrected  for  tunnel  effects,  (Cland  Cd) 

^i?kcicickc*icic*kcicickckcicic*kcicickc*icic*kcicickc*icickckcicickc*icic*kcicickcjcicic*kcic'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k-k'k'k'k'k 


C  L  u  =  Forces 

w (3,  :  ) 

./ 

(q  Corrected' 

. *  Wing  Area) 

%Keuthe  &  Chow  pg  178 

C  D  u  =  Forces 

w (1,  :  ) 

./ 

(q  Corrected' 

. *  Wing  Area) 

C  Y  u  =  Forces 

w (2,  :  ) 

./ 

(q  Corrected' 

. *  Wing  Area) 

Coefficients  = 

[C  L  u; 

C 

Du;  C  Y  u]  '  ; 

%  Ave_Cl  =  mean (Coefficients (:, 1 )) ; 

%  Ave_Cd  =  mean  (Coefficients (:, 2) ) ; 

end 

%******************  Drag  Coefficient  Correction  *************************** 

SJj' k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 


C_D_o  =  min (Coefficients (:, 2 )) ; 

C_L_u_sqrd  =  Coef ficients ( : , 1 )  . A2  ; 

Delta_C_D_w  =  ((delta  *  Wing_Area)  /  X_Section)  .*  C_L_u_sqrd; 

C_D_Corrected  =  C_D_u '  +  Delta_C_D_w; 

~k 
~k 
' k 


S^'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k 

%*************  Angle  of  Attack  due  to  upwash  Correction  ****************** 


alpha  =  sample_data ( : , 1 ) ; 

Delta_alpha_w  =  ((delta  *  Wing_Area)  /  X_Section)  .*  (57.3  *  C_L_u) ; 
alpha_Corrected  =  alpha  +  Delta_alpha_w ' ; 

^i?kcicickcj?icic*kcicickckcicicjckcicickcjcicic*kcicickckcicic*kcicickcj?icic*kcicickcjcicifjcjci<:ickcj?icic*kc'k'k'k'k'k'k'k'k'k'k'k-k'k'k'k'k 

%*************  pitching  Moment  Correction  ********************************* 


Cl_ 

w 

.eg  = 

Moments 

w 

_cg_ 

u (1,  : ) 

./ 

(q 

Corrected ' 

.  *  (Wing 

Area  * 

Span*12) 

Cm 

w 

eg  u 

=  Moments 

w 

_cg_ 

u (2 ,  : ) 

./ 

(q 

Corrected ' 

.  *  (Wing 

Area  * 

c  bar*12 

Cn_ 

w 

.eg  = 

Moments 

w 

_cg_ 

u (3,  : ) 

./ 

(q 

Corrected ' 

.  *  (Wing 

Area  * 

Span*12) 

Cm_w_cg_corrected  =  Cm_w_cg_u;  %-  Delta_C_m_cg_t ' ; 

Corrected_Moment_Coef ficients  =  [Cl_w_cg'  Cm_w_cg_corrected '  Cn_w_cg']; 
%**Obtaining  moments  coefficients  corrected  about  the  center  of  the  balance 

S^'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

Cl_w_bc  =  Moments_w_bc ( 1 , :)  ./  (q_Corrected '  .*  (Wing_Area  *  Span  *  12)); 

Cm_w_bc_u  =  Moments_w_bc ( 2 , : )  ./  (q_Corrected '  .  *  (Wing_Area  *  c_bar  *  12)); 

Cn_w_bc  =  Moments_w_bc (3, :)  ./  (q_Corrected '  .*  (Wing_Area  *  Span  *  12)); 

Cm_w_bc_corrected=Cm_w_bc_u; 

Corrected_Moment_Coef ficients_bc= [Cl_w_bc '  Cm_w_bc_corrected '  Cn_w_bc ' ] ; 

^•k  k  ~k  ★  •*•  •*•  ★  ★  ~k  k  ~k  ~k  k  k  ★  ★  •*•  *  ★  ★  *  *  ★  ★  *  *  ★  ★  •*•  *  ★  ★  •*•  •*•  ★  ★  •*•  *  ★  ★  ~k  k  ~k  ~k  k  k  ~k  ~k  k  k  ~k  ~k  k  k  ★  ★  •*•  *  ★  ★  •*•  -*•  ★  ★  *  ~k  ~k  ~k  k  k  ~k  ~k  k  k 

%****************** Output  Variables  Formatting***************************** 
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Q''k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 


alpha  =  sample_data ( : , 1 ) ; 

C_Y\r ' ) ; 

YY= [Flight_Parameters  (Wind_Speed_Corrected  .*  (3600/5280))  alpha_Corrected 
C_L_u '  C_D_Corrected  Corrected_Moment_Coef f icients  C_Y_u ' ] %pressure ] 
wklwr ite ( ' output . xls '  ,  YY,  2, 1) 

Max_Cl  =  max (Coefficients (:,1)); 


figure; 

plot (alpha_Corrected, Coefficients ( : , 1) , 'b . - . ' ) ; 
legend ('100  mph  ' , ' \beta  \approx  0\circ'); 
grid  on; 

title  (  ' \it  C_L  vs  \alpha  no  winglets  (50 

mph ) 1 , ' FontWeight ' , ' bold ' , ' FontSize ' , 11 ) ;  xlabel (' Angle  of  Attack  (\alpha)') 
ylabel ( 'Lift  Coefficient  (CL) ' ) ; 
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In  an  effort  to  increase  the  range  of  missiles  and  guided  bombs,  the  USAF  is  looking  at  options  for  new  wing  configurations. 
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less  drag  for  a  given  lift  at  both  supersonic  cruise  and  subsonic  loiter  conditions.  This  experiment  focused  on  the  low  speed 
performance  of  a  missile  model  with  an  oblique  wing.  The  wing  was  tested  at  seven  different  sweep  angles  and  at  two 
different  speeds.  In  order  to  simulate  the  missile  dropping  from  an  aircraft  the  model  was  inverted  over  a  stationary  ground 
plane  in  the  tunnel  and  tested  at  the  same  wing  obliquity  angles.  Stalling  was  found  at  certain  conditions  including  sweep 
angles  of  0,  15,  and  30  deg.  The  ground  plane  was  shown  to  result  in  an  increase  in  lift  as  well  as  an  increase  in  drag.  The 
ground  plane  was  also  shown  to  add  more  longitudinal  stability,  thus  making  the  missile  better  performing  when  dropped  from 
an  aircraft. 
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